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Electronic Structure of transition metals Fe, Ni and Cu in the GW approximation 
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The quasiparticle band structures of 3d transition metals, ferromagnetic Fe, Ni and param- 
agnetic Cu, are calculated by the GW approximation. The width of occupied 3d valence band, 
which is overestimated in the LSDA, is in good agreement with experimental observation. How- 
ever the exchange splitting and satellite in spectra are not reproduced and it is required to go 
beyond the GW approximation. The effects of static screening and dynamical correlation are 
discussed in detail in comparison with the results of the static COHSEX approximation. The 
dynamical screening effects are important for band width narrowing. 

KEYWORDS: Density functional theory, GW approximation, LMTO, product-basis, transition metal, dynamical 
screening 



%4 , 1. Introduction 



Lattice structure, lattice constants and bulk moduli 
in 3c? transition metal are well described by the local- 
spin-density approximation (LSDA) 1 ' 2 ) or the general- 
ized gradient approximation (GGA). 3,4 ) However, the 
occupied 3c? band width is too broad, and the exchange 
splitting is overestimated. The good agreement is essen- 
tially related with the property of the ground-state, and 
the discrepancies are associated with excitation proper- 
ties. 

The GW approximation (GWA) is based on the many- 
body perturbation theory 5-7 ) and can describe the quasi- 
particle property. The self-energy of GWA is the first 
term in a series expansion of dynamical correlation and 
it is treated by the random-phase approximation (RPA) . 

The plane wave basis set based on the pseudopoten- 
tial method is used in many GW calculations. In sim- 
ple metals and semiconductors, the single plasmon peak 
is often assumed within the plane wave framework (the 
plasmon pole approximation). 8 ) However the plasmon 
peak of transition metal cannot be well-defined isolated 
peak due to interband transition in the same energy re- 
gion. The transition metal has strong atomic potential 
for 3d electrons, the 3d orbital is localized and the plane 
wave formalism cannot be applied. Moreover it is essen- 
tially important to include core electrons in many cases. 
Therefore the plasmon peak approximation is not appli- 
cable to the dielectric function of transition metals and 
all-electron calculation and localized orbital basis set are 
needed. 

In this paper, the GW method based on the linear 
muffin-tin orbital (LMTO) method 9 ) and the product- 
basis method 10 ) are applied to the series of transition 
metal. There is a numerical difficulty in the k-point sum- 
mation of self-energy with the momentum transfer q = 0. 
This summation is treated by the offset method, 11 ) and 
test calculation of the exchange energy in the electron 
gas is performed. The paper is organized as follows. The 
theoretical framework is described in § 2. The numeri- 
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cal technique and test calculation of electron gas are also 
given in this section. The results for these systems and 
detailed discussion are presented in § 3. Finally, in § 4 
we present our summary 

2. Theory 

2.1 GW approximation 

In the GWA the self-energy is replaced by the lowest 
order term of the expansion as E(l, 2) = iG(l, 2)W(1, 2). 
G is the one particle Green function and the dynamically 
screened interaction W is defined by 



W(l,2)= / d(3)e- 1 (l,3) U (3,2) 



(1) 



«(1,2)+ / d(34Ml,3) X °(3,4)W/(4,2), (2) 



where e 1 is the inverse dynamical dielectric function, v 
is the bare Coulomb potential and x° is the irreducible 
polarization function x°(l,2) = -iG(l, 2)G(2, 1). Here 
we use an abbreviated notation (1) = (ri,<7i,ti) and 
v(l,2) — v(ri, r 2 )S(ti — t 2 ). Equation (2) is treated by 
the RPA. 

We adopt the LSDA Hamiltonian to be the unper- 
turbed one H° = T + V H + V£§ DA . Here T is the kinetic 
energy, V H is the Hartree potential, and V£$ DA is the 
exchange-correlation potential in the LSDA. We presume 
the wavefunctions {V'kn(r)} of the LSDA to be a reason- 
ably good starting wavefunctions. Then the self-energy 
can be written by three terms as AS = E K +E C — V£g DA , 
where E x (= iGv) is the exchange part (the Fock term) 
and £ c (= iGW c ) is the dynamical correlation part. W° 
is the second term in eq.(2). The quasiparticle energy is 
given as 



(3) 



where £kn is the LSDA eigenvalue. The self-energy is 
A£ kn (e k „) = (i/>k n \Z x + S c (e k „) - V£§ DA \t/tkn) and the 
renormalization factor is Z^ n = (1 — 9A£k n (ekn)/9w) . 
The renormalization factor Z^ n is a measure of the oc- 
cupation number and should equal to the discontinuity 
of occupation number at the Fermi energy. Therefore 
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Fig. 1. Exchange energy S(y) of the electron g function of 

y = k/kp. The solid lines for the exact result, the open triangles 
for the simple summation and the closed circles for the offset 
method. The numbers of the mesh points in the Brillouin zone 
are (a) 64 = (4 X 4 X 4) and (b) 512 = (8 X 8 X 8). 



it should satisfy the condition Z^ n < 1. In the present 
work we perform one iteration calculation without self- 
consistency. 12 ) 

2.2 LMTO minimal basis set and product-basis 
Because the plane wave basis becomes very costly for 

systems containing 3d electrons, the LMTO method 9 ^ is 
more appropriate. We use the LMTO basis set XRLv(r) 
within the atomic sphere approximation (ASA) for the 
LSDA calculation. Here L is angler momentum L = 
(I, m). The LMTO can be expanded by the muffin-tin 
orbital 0RLi/(r) and it's energy derivative 4>rlu(^)- 
The functional space of basis for S is spanned as 

{£} = {^} = {xx} = {H>} + {H>} + {#} ■ (4) 

In fact the mixing coefficients of <p to <j> are less than 0.1 
for the most part and the norm of </>, {4>4>), is 0.1 ~ 0.3 
even in the largest case, so the terms including </> can 
be dropped out. More detailed description is shown in 
rcf. 13. 

2.3 Numerical technique 

The Coulomb matrix v(q) has a singularity at q = as 
F(q) = l/|q| 2 - The integration of v(q) over the Brillouin 
zone does not diverge but special cares are needed not 
only for the q = term but for small finite q. For 
a choice of the discrete points near q = 0, we use the 
offset T-point method, n ) where the integration of F(q) 
over the Brillouin zone can be performed analytically and 
the offsetted points Q's are chosen near q = so as to 
satisfy a relation 



/ 

JB.Z. 



F(q)dq 



Q 



(5) 



Here k's are the discrete mesh points in the Brillouin 
zone. 

The exchange energy of the electron gas system is 
given as a function of a wave vector k as Y. x (k) = 



^S(y), where S(y) = - (l + In | i±* |) , y = £ 
and kp is the Fermi wave vector. The empty lattice cal- 
culation is done with spdf orbitals in the LMTO method. 
We calculate the exchange energy of the electron gas in 
a fee lattice with a lattice constant a = 6.824ao which 
corresponds to the fee copper and do is the Bohr radius. 
The corresponding electron gas parameter is r s = 2.6668. 
The calculated S(y), by the simple summation, by the 
offset method and by the exact S(y), are shown in Fig. 1. 
In the simple summation, the diverging term l/|q| 2 is 
simply averaged inside a sphere of a volume equal to 
that of one k-mesh point. The number of k-mesh of 
the Brillouin-zone in the calculation is (a) 64 points 
(4x4x4) and (b) 512 points (8x8x8). In case of 
512 points of these structure, eight offsetted points Q 
arc ^(±0.038, ±0.038, ±0.038). The derivative of S(y) 
has a logarithmic singularity at the Fermi energy (y = 1). 
Unphysical gap still remains at the Fermi energy in the 
simple summation of these examples. A large number 
of k points is necessary for a convergence in the sim- 
ple summation. But the offset method can reduce the 
number of k point for rapid convergence even in case of 
small number of mesh points. The careful treatment of 
the Coulomb matrix at or near q = is very crucial near 
the band gap or the Fermi energy. 

3. Results and Discussions 

In the calculation of LSDA, the lattice structure and 
constants of Fe, Ni and Cu are bee and a = 2.87 A, fee 
and 3.52 A, fee and 3.61 A, respectively. 14 ^ The band 
structures of Fe, Ni and Cu, calculated both in the LSDA 
and the GWA, are shown in Fig. 2 along high symmetric 
lines. The localized 3d orbital has a weak hybridization 
with the extended 4s, Ap orbitals and is below Fermi 
energy. 

The magnetic moment, the exchange splitting and the 
band width of the occupied id valence bands in the GWA 
are summarized in Table I, in comparison with those 
by the LSDA and the static COHSEX approximation. 6 ) 
Our results of Ni are in good agreement with those of 
the previous GW calculation. 15 ) The spectral function 
A(uj) = -(l/7r)ImTrG(w) is shown in Fig. 3. 



Table I. The Magnetic moment fJ. sp i n (fJ-B) °f Fc an d Ni, the ex- 
change splitting SEex (eV) of Fe and Ni, the band width of the 
occupied 3d valence bands Wd,occ ( e V) of Fe, Ni and Cu. 
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Fig. 2. The energy bands of Fe, Ni and Cu, calculated by the LSDA (dotted lines) and the GWA (solid lines) along high symmetric 
lines. The closed circles are the calculated points in the GWA. In Fe and Ni, left side and right side show minority spin and majority 
spin, respectively. The high symmetric points arc V = (0, 0, 0), H= (1, 0, 0), N= (1/2, 1/2, 0) and P= (1/2, 1/2, 1/2) in bec lattice (Fc), 
and L= (1/2, 1/2, 1/2), T = (0, 0, 0) and X= (1, 0, 0) in fee lattice (Ni and Cu). Fermi energy is set to zero. (E F =0) 



The occupied 3d valence band width of transition met- 
als Fe, Ni and Cu in the LSDA is overestimated in com- 
parison with experimental observation, especially in Ni. 
The valence band width is in reasonably agreement with 
experiment in the GWA. In Fig. 3, the band narrowing 
occurs in the occupied valence band of both the major- 
ity and the minority spin in Fe and Ni. But the width 
of unoccupied 3d band of Fc is unchanged in the GWA. 
The source of band narrowing is the screening for the 
valence electrons. In the spectral function of the GWA, 
the plasmon-like excitation appears around 30 eV above 
and below the Fermi energy and also the long tail ex- 
tends over wide lower energy region. The intensity of 
spectrum is totally suppressed by the excitations in wide 
energy region. The intensity of the GW spectrum is ac- 
tually reduced by a factor of Z kn and, in Fig. 3, the 
reduction factor is Zd « 0.5-0.6. In the Hartree Fock 
(HF) approximation which includes no screening effects, 
the band width is overestimated. No screening in the 
HF gives zero density of states at the Fermi level in the 
electron gas and also gives overestimated band gap in 
insulators and semiconductors. In the static COHSEX 
approximation, which includes static screening, the band 
width is much smaller than the HF results, and is almost 
the same as the one in the LSDA or still wider. Moreover 
the 4s state is located too much deep because it exists far 
from Fermi level. The band width in the GWA is in good 
agreement with experiment. We can see that the dynam- 
ical correlation effect is important for the band width in 
the transition metals from the comparison between the 
GWA and the static COHSEX approximation. 

The magnetic moment (i sp i n of Fe and Ni is almost 
the same as the result of LSDA and is in good agree- 
ment with experiment. The difference of the exchange 
splitting SE ex between the LSDA and experiments in Fe 
(~ 30%) is smaller than that in Ni (~ 50%). SE ex of 
Fe becomes close to the experimental value in the GWA. 
In Ni, the discrepancy of SE ex is not improved by the 
GWA. In the HF, 5E ex is overestimated. The screening 
effects of correlation term S c in the GWA or the static 
COHSEX approximation reduce SE ex of HF. However 



the GWA only includes long-range correlation effects, 
and cannot describe short-range effects such as electron- 
electron or hole-hole scattering process. Higher order 
diagrams (e.g. vertex corrections) is needed for electron- 
electron and hole-hole scattering. Especially two-hole 
bound states are very important to the exchange split- 
ting and the satellite structure of spectrum if on-site 
Coulomb interaction between d electrons is large. 27, 2S ) 
The effective Coulomb interaction is obtained from an 
analysis of Auger spectra, Ni is U « 4.0 eV and Fe is 
U w 1.0 eV. 29 ) The discrepancy between experiments 
and the GWA in Ni is caused by the short-range correla- 
tion effects from large U. The clear satellite in Ni can be 
also explained, and the GWA cannot reproduce it. Since 
the 3d band is full in Cu, there is no hole-hole correlation 
and the GWA can work quite well. 

The rcnormalization factor of transition metal 3d 
states is Z = 0.52-0.58 in Fc, Z = 0.48-0.53 in Ni and 
Z = 0.53-0.66 in Cu. Z of 4s states is about 0.7-0.8 
in these systems. Those results of the renormalization 
factor show that the interaction between 3d electrons is 
large, and the correlation in Ni is strongest, which is 
consistent with the large Coulomb interaction U. 

The static screened d-d Coulomb interaction (4>d4>d\ 
W(uj = 0)|(/>d</>(j) is calculated to be about 1.4 eV, 1.2 eV 
and 3.9 eV in Fe, Ni and Cu, respectively. The bare 
Coulomb interaction {4>d4>d\v\4>d4>d) is actually 23.7 eV, 
25.9 eV, and 27.4 eV in Fc, Ni and Cu. This values 
become larger with increasing the number of 3d occupa- 
tion. Then the correlation term ((j)d<f>d\W c (w = Q)\(j)d<i>d) 
is -22.3 eV, -24.7 eV and -23.5 eV in Fc, Ni and Cu. 
Therefore, the correlation effects and the static screening 
are quite important in transition metals. The screened 
correlation of Ni is largest and this is consistent with 
the smallest Z in these systems. We should mention 
that the d-d Coulomb interaction (<fid ( j ) d\W(0)\(j>d(j>d} is 
different from the Hubbard U evaluated from the con- 
strained LSDA, which includes only screening by on-site 
d electrons. The term (4>d4>d\W{0)\4>d4>d) includes the 
screening effects by both on-site and off-site electrons. 30 - ) 

In the transition metals, the Hubbard U parameter is 
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Fig. 3. The spectral function A(lu) in the transition metals Fe, 
Ni and Cu by the GWA (bottom), the LSDA (middle) and the 
static COHSEX approximation (top). In Fe and Ni, the solid 
lines show majority spin and the dotted lines show minority spin, 
respectively. 



overestimated within the constrained LSDA, for example 
U ~ 6 eV for Fe, due to incomplete metallic screening 
in the LSDA. 31 ) We should mention that our value of Ni 
with offset method is smaller than the previous estimate 
{{(t>d<t>d\W{<S)\<t>d<t>d) =2.2eV). 3 °) The discrepancy maybe 
caused by the absence of the present offset method since 
we also obtained the value without the offset method 
similarly to be previous one. 

4. Summary 

In this paper the GW approximation is applied to fer- 
romagnetic transition metals Fe and Ni, and paramag- 
netic Cu. We showed that the occupied 3d band width of 
transition metal is improved within our GW calculation. 

We also investigated the effects of dynamical screening 
by comparison between the GWA and the static COH- 
SEX approximation and showed the crucial role of the 
dynamical correlation for band width. 

The self-energy is discussed systematically. The renor- 
malization factor Z showed that the interaction between 
d-electrons in Ni is larger than that in Fe, and this is con- 
sistent with the interaction strength from Auger spectra. 
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